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Abstract 


During the 1998 and 2001 hurricane seasons of the western Atlantic Ocean and 
Gulf of Mexico, the Advanced Microwave Precipitation Radiometer (AMPR), the ER-2 
Doppler (EDOP) radar, and the Lightning Instrument Package (LIP) were flown aboard 
the National Aeronautics and Space Administration ER-2 high altitude aircraft as part of 
the Third Convection and Moisture Experiment (CAMEX-3) and the Fourth Convection 
and Moisture Experiment (CAMEX-4). Several hurricanes, tropical storms, and other 
precipitation systems were sampled during these experiments. An oceanic rainfall 
screening technique has been developed using AMPR passive microwave observations of 
these systems collected at frequencies of 10.7, 19.35, 37.1, and 85.5 GHz. This 
technique combines the information content of the four AMPR frequencies regarding the 
gross vertical structure of hydrometeors into an intuitive and easily executable 
precipitation mapping format. The results have been verified using vertical profiles of 
EDOP reflectivity and lower altitude horizontal reflectivity scans collected by the 
National Oceanic and Atmospheric Administration WP-3D Orion radar. Matching the 
rainfall classification results with coincident electric field information collected by the LIP 
readily identifies convective rain regions within the precipitation fields. This technique 
shows promise as a real-time research and analysis tool for monitoring vertical updraft 
strength and convective intensity from airborne platforms such as remotely operated or 
uninhabited aerial vehicles. The technique is analyzed and discussed for a wide variety of 
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precipitation types using the 26 August 1998 observations of Hurricane Bonnie near 
landfall. 
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1. Introduction 


A crucial need exists to understand and map the precipitation types, patterns, and 
variations of a tropical cyclone (TC) in order to develop better skill in quantitative 
precipitation estimation necessary for more accurate forecasts of rainfall impacts during 
landfalling tropical storms and hurricanes. Identification of rainfall patterns, vertical 
hydrometeor profiles, and corresponding vertical motions are also necessary for defining 
latent heat profiles and regions of convective strength, which in turn can be used to 
improve hurricane intensity change forecasting as well as general numerical weather 
prediction. 

Efforts to identify precipitation characteristics of tropical oceanic convective 
systems, in general, and tropical cyclones, in particular, have included extensive analysis 
of microwave remote sensing information. Early works included Rodgers and Adler 
(1981) who used passive microwave Nimbus 5 information to document rain patterns of 
eastern and western Pacific tropical cyclones, and Wilheit et al. (1982) who discussed the 
correlation of increasing 19.35 GHz and decreasing 92 GHz passive microwave brightness 
temperatures to increasing rain rates using aircraft observations of Tropical Storm Cora. 
Jorgensen (1984) used aircraft radar reflectivities to identify meso scale and convective- 
scale features of mature hurricanes while Marks (1985) used aircraft radar reflectivities 
collected during Hurricane Allen in 1980 to examine the relationship of storm intensity 
change to rain rate and total rainfall. Burpee and Black (1989) used National Weather 
Service WSR-57 radar reflectivities of Hurricanes Alicia in 1983 and Elena in 1985 to 


5 


identify asymmetric trends of precipitation distribution in hurricane eyewall and rainband 
regions. 

More recent work include Spencer et al. (1994) and McGaughey et al. (1996) who 
presented different aspects of passive microwave aircraft observations of Tropical 
Cyclone Oliver collected by the Advanced Microwave Precipitation Radiometer (AMPR) 
during 4-9 February 1993 in the western Pacific as part of the Tropical Ocean and Global 
Atmosphere Coupled Ocean - Atmosphere Response Experiment (TOGA-COARE). 
Spencer et al. (1994) presented the first high-resolution, multi-frequency passive 
microwave imagery of a TC during their description of the AMPR sampling capabilities. 
McGaughey et al. (1996) explored multi-frequency passive microwave signatures of 
tropical oceanic precipitation systems. They explained the spatial shift of lower altitude 
rain emission microwave signatures from higher altitude ice scattering microwave 
signatures as a result of the tilt with height of convective elements in the eyewall of TC 
Oliver. Tilted convective elements within Hurricane Bonnie on 25 August 1998 in the 
western Atlantic Ocean have also been identified by Hong et al. (2000) using Tropical 
Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) and Special Sensor 
Microwave/Imager (SSM/I) satellite observations. Cecil and Zipser (1999) examined the 
relationship of satellite observations of passive 85.5 GHz ice scattering signatures and 
lightning in TC eyewalls and rainbands to future TC intensity change. Cecil and Zipser 
(2002) examined relationships between satellite passive microwave, radar, lightning and 
inferred microphysical characteristics of eyewalls and rainbands. Simpson et al. (1998) 
explored cloud electrification and lightning linked to the vertical radar structure and other 
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features of the clouds in a study of cyclogenesis in TOGA-COARE associated with TC 
Oliver. 

In general, satellite, aircraft and ground-based information have their own 
advantages and disadvantages for TC study. Aircraft and ground-based instrumentation 
provide detailed information of TC features but may be limited in complete spatial 
coverage of a TC. Conversely, satellites are able to provide more complete spatial 
coverage of a TC in a wide variety of global locations including remote areas inconvenient 
for aircraft study. However, the coarser spatial resolution of some satellite observations 
may not be able to discern convective-scale features. Temporal sampling from low earth 
orbit is limited to every three hours at best. These lower spatial and temporal 
resolutions from satellites are disadvantages for understanding TC structure and intensity 
change. Additionally, many satellite algorithms still require aircraft or ground-based 
information for validation purposes. 

In an effort to address a need for detailed information regarding the variability of 
TC characteristics within the global water cycle and to validate satellite moisture 
measurements of the TRMM, a comprehensive volume of information using spacebome, 
airborne, and ground-based instrumentation was collected during the Third Convection 
and Moisture Experiment (C AMEX-3) in 1998 and the Fourth Convection and Moisture 
Experiment (C AMEX-4) in 2001. Several hurricanes, tropical storms, and other 
precipitation systems in the western Atlantic and Gulf of Mexico were sampled during 
these field campaigns. Both experiments were sponsored by the Earth Science Enterprise 
of the National Aeronautics and Space Administration (NASA). A collaborative 
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partnership with the Hurricane Research Division of the National Oceanic and 
Atmospheric Administration (NOAA) provided opportunities for joint missions during 
C AMEX-3 and CAMEX-4 using NOAA WP-3D Orion (P-3) aircraft and 
instrumentation. 

This study demonstrates a methodology to merge the information content of 
several CAMEX remote sensors into a format that highlights the type and convective 
strength of the TC precipitation elements sampled. In particular, this study examines 
information collected by the AMPR, the ER-2 Doppler (EDOP) radar, and the Lightning 
Instrument Package (LIP), which were deployed aboard the NASA ER-2 high altitude 
aircraft. Descriptions of the instruments are presented in Section 2. A rainfall screening 
technique has been developed using AMPR passive microwave observations of these 
tropical cyclones and other precipitation systems sampled during other field 
opportunities. The rainfall classification was then verified using vertical profiles of 
EDOP reflectivity and lower altitude horizontal reflectivity scans collected by the NOAA 
P-3 radar. Matching the rainfall classification results with coincident electric field 
information collected by the LIP readily identifies convective rain regions within TC 
precipitation fields. The AMPR screening method as well as the EDOP verification and 
coincident LIP data are presented in Section 3. An illustration is presented in Section 4 
using a portion of Hurricane Bonnie data collected on 26 August 1998 that displays a 
wide variety of precipitation structures. This technique shows promise as a real-time 
analysis tool for monitoring vertical updraft strength and convective intensity from a 
remotely operated or uninhabited aerial vehicle. These types of vehicles are likely 
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components of a future network of spacebome, suborbital, and ground-based Earth 
observing platforms combining the advantages of each platform for flexible, adaptive 
sampling of critical weather events such as tropical storms and hurricanes. 

2. Instrument Descriptions 

a. The Advanced Microwave Precipitation Radiometer (AMPR) 

The AMPR is a total power passive microwave radiometer producing calibrated 
brightness temperatures (TB) at 10.7, 19.35, 37.1, and 85.5 GHz. These frequencies are 
sensitive to the emission and scattering of precipitation-size ice, liquid water, and water 
vapor. The AMPR performs a 90° cross-track data scan perpendicular to the direction of 
aircraft motion, with full vertical polarization at -45° and frill horizontal polarization at 
+45°. The polarization across the scan is mixed as a function of sin 2 , giving an equal V-H 
mixture at 0° (aircraft nadir). A frill calibration is made every fifth scan using hot and cold 
blackbodies. From a typical ER-2 flight altitude of ~20 km, surface footprint sizes range 
from 640 m (85.5 GHz) to 2.8 km (10.7 GHz). A more complete description of the 
instrument may be found in Spencer et al. (1994). 

An example of AMPR TB imagery is presented in Fig. 1 for a portion of Hurricane 
Bonnie on 26 August 1 998. The false color scale chosen for this imagery ranges from 
magenta and blue colors for cold TB to yellow and red colors for warm TB. In general, 
land surfaces and rain are radiometrically warm while the ocean is radiometrically cold due 
to their respective microwave emission properties. Precipitation-ice may also appear as a 
radiometrically cold signature for a given frequency because ice tends to scatter upwelling 
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microwave energy out of the instrument field of view unless the diameter of ice particles 
is small compared to the wavelength. See Wilheit et al. (1977), Wu and Weinman (1984), 
Wilheit (1986) and Spencer et al. (1989) for further explanation. 

In Fig. 1, the 10.7 GHz information delineate rain and rain-free regions with wanner 
TB representing increasing rain rates. The 19.35 GHz information also serve this 
purpose but with a smaller dynamic range than the 10.7 GHz. (Smith et al. (1994). The 
19.35 GHz frequency is more sensitive to clouds and ice than the 10.7 GHz. In die 
image, the 37.1 GHz information displays a similar sensitivity to clouds as the 19.35 
GHz frequency, but the 37.1 GHz frequency is more sensitive to smaller size ice. The 
coldest 37.1 GHz TB in this image represent clear sky regions over a radiometrically cold 
ocean background such as the eye near 33.2N and 77.8W. In the 85.5 GHz image, water 
vapor, clouds, and smaller ice are very noticeable. Precipitation ice is colder than the 
surrounding rain and cloud. Regions of ‘blues’ and ‘greens’ (e.g., at 32.6N and 77.4W) 
indicate large quantity and/or mass of ice. Note that the eye is almost obscured at this 
frequency probably due to thin clouds and/or high water vapor content. The spray of 
small dots found in Fig. 1 corresponds to geolocation during aircraft turns. 

b. The ER-2 Doppler Radar (EDOP) 

The EDOP, operating at 9.6 GHz, provides high resolution (i.e., 37.5 m vertical, with 
the horizontal footprint varying from about 400 m at the tropopause level to 1.2 km at 
the surface) time-height sections of reflectivity and vertical hydrometeor velocity in the 
vertical plane mapped out by the ER-2 (Heymsfield et al., 2001). Vertical air motions 


1ft 


are also retrieved when the hydrometeor fall speed and aircraft motions are removed. 
EDOP reflectivities are calibrated to within approximately 1 dBZ. These have been 
verified using various approaches including the use of the ocean backscatter and 
comparison with the TRMM Precipitation Radar and ground-based radars. The minimal 
detectable signal is approximately -15 dBZ near cirrus cloud top and 0 dBZ near the 
surface. A more complete description of the EDOP and its other capabilities (e.g., ability 
to measure linear depolarization ratio, dual Doppler retrieval along the flight track) may 
be found in Heymsfield et al. (1996). 

c. The Lightning Instrument Package (LIP) 

The LIP consists of eight state-of-the-art, low-noise, high dynamic range electric field 
mills on the aircraft (three mills per instrument superpod mounted on each wing and two 
on the fuselage). With these sensors, the full vector components of the atmospheric 
electric field (i.e., E*, Ey, E z ) are obtained, providing detailed information about the electric 
structure within and around the storms overflown. The field mills measure the 
components of the electric field over a wide dynamic range extending from fair weather 
electric fields (i.e., a few V m' 1 ) to large thunderstorm fields (i.e., tens of kV m' 1 ). The set 
of equations that relate the field mill outputs to the atmospheric electric field is 
determined by an iterative calibration process (Mach and Koshak 2003). Total lightning 
(i.e., intracloud and cloud-to-ground) is identified from the abrupt changes in the electric 
field data. 

3. Method of analysb 
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a Hydrometeor classification using AMP R 

An AMPR Precipitation Index (API) has been developed that utilizes the brightness 
temperature information for precipitation and clouds over an ocean background at the four 
frequencies to produce a single index value at each AMPR footprint. The goal of the API 
development is to combine the information content of the four frequencies into an 
intuitive format that readily identifies the gross vertical structure of the hydrometeors at a 
given pixel location. The method is not meant to be a replacement for the iterative 
hydrometeor retrieval of Skofronick-Jackson et al. (2002) or the texture-polarization 
method of Olson et al. (2001). The Skofronick-Jackson retrieval method uses nadir 
aircraft observations of active and passive microwave sensors and a cloud resolving model 
to deduce vertical content and particle size distribution. The Olson method is a satellite 
technique for conically scanning passive microwave radiometers that estimates the area 
coverage of convective and stratiform precipitation using 85.5 GHz polarization 
information and lower frequency texture data correlating local maximum signatures to 
neighboring footprints. The API, on the other hand, is presented as an alternate approach 
that is a computationally easy method to map and identify precipitation type for the 
extensive precipitation data sets collected by the AMPR during CAMEX-3 and 
CAMEX-4 or future data sets requiring real-time analysis. The API technique presented 
here is dependent upon the scanning strategy of the AMPR, but could be readily adapted 
for other radiometer scanning strategies or expanded to include information from other 


In general, the API reflects the magnitude (mass) of liquid water and precipitation- 
sized ice aloft It is based on physical concepts of microwave rain emission and ice 
scattering discussed in the literature (e.g., Wilheit et al. 1977, Wu and Weinman 1984, 
Wilheit 1986, Spencer et al. 1989, Smith et al. (1994)). The indices are listed in Table 1 
along with color codes (for imagery), descriptors, and estimated rain rates. These rain 
rates are presented for illustrative purposes only. The conversion of TB to quantitative 
rain rate is beyond the scope of this paper. Instead, the emphasis is on initial 
precipitation screening following the example of Ferraro et al. (1998). 

The API is designed for ocean-only cases that have been screened for large aircraft 
pitch, roll, and altitude variations that greatly influence the TB values. All pixels within 
3.2 km of land are not used to avoid contamination horn varying land surface emissions. 
An API value of 0 indicates that no clouds or rain are detected by AMPR. API values of 
1 or 2 suggest the presence of clouds or veiy light rain, with microwave emissions 
exceeding 190 K (19.35 GHz) or 260 K (85.5 GHz). The AMPR frequencies are not well 
suited to detecting some cloud types; so we do not interpret an API value of 0 as 
indicating truly clear skies. 

For API rain values of 3 through 1 8, six TB rain emission tests and four TB ice 
scattering tests are performed. First a rain / no-rain emission test is performed based on 
10.7 GHz and 37.1 GHz TB. The 10.7 GHz channel is used because it is most sensitive 
to emission by liquid rain and least sensitive to scattering by ice. The 37.1 GHz channel 
is used to resolve small-scale features due to its higher spatial resolution. The thresholds 
for this test are TB10 > 160 K or TB37 > 215 K at nadir. They vary across the scan by 
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up to 53 K for the 10.7 GHz test and 30 K for the 37.1 GHz test to account for the 
AMPR rotating polarization. After this test, the remaining five emission tests check for 
TB10 exceeding 175, 200, 225, 250, 275 K, indicating increasing liquid water mass and, 
thus, inferring increasing rain rates. 

The scattering tests indicate which AMPR wavelengths are being scattered by ice. 
In general, the larger the ice particles that are present, die longer the wavelength that will 
be scattered. This relationship can be used as a surrogate indicator of vigorous 
convection. For liquid rainfall in the absence of appreciable precipitation-sized ice, die 
higher frequency channels usually have greater TB than the lower frequency channels. In 
such cases, the precipitation index has values of 3 through 5. These are shown in blue in 
Table 1, with darker blues indicating greater 10.7 GHz emission (i.e., more rain). 
However, if TB85 is less than TB37 and less than a threshold of 275 K, it is likely 
displaying the effects of ice scattering. If only the 85.5 GHz channel is scattered (i.e., ice 
is large enough to scatter 3.5 mm wavelength radiation), an index of 6 through 1 0 is 
assigned (green shades). If TB37 is also less than TB19, we interpret this to mean the ice 
is large enough to scatter the 37.1 GHz channel (8.1 mm wavelength), and an index of 1 1 
to 1 5 is assigned (shades of yellow/red). If the ice present is large enough to scatter the 
19.35 GHz channel (1.6 cm wavelength), then TB19 is decreased to a value less than 
TB10. An index of 16 to 18 is assigned (shades of violet/purple). The strongest 
convection (with the largest ice) also scatters the 10.7 GHz channel, but this is not 
incorporated into the algorithm because a lower frequency would be required in the 
current framework. 
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b. Verification using EDOP information 

The API has been compared with EDOP reflectivity profiles for the tropical cyclone 
and precipitation cases collected during CAMEX-3 and CAMEX-4. This includes 
missions over Hurricane Bonnie (23, 26 August 1998), Hurricane Earl (2 September 
1998), Hurricane Georges (21, 22, 25, 27 September 1998), Hurricane Erin (10 September 
2001), Hurricane Humberto (22-24 September 2001), and other convective systems near 
Florida (5, 17 September 1998; 9, 19 September 2001). For each AMPR scan, the 
precipitation indices of the middle two pixels (i.e., those nearest nadir) are matched with 
the simultaneous nadir reflectivity profile from EDOP. This yields -80,000 realizations 
of the vertical profiles associated with the AMPR Precipitation Indices. 

From the set of observed reflectivity profiles, a characteristic (median) profile is 
assigned to each API. These characteristic profiles are shown in Fig.2, and are used to 
produce simulated radar reflectivity from the AMPR measurements. The variability 
about each characteristic profile is assessed using cumulative density function (CDF) of 
reflectivity. Such CDFs for some of the most common API values are shown in the 
panels of Fig . 3. Although beyond the scope of this paper, the precipitation index can 
be applied to various problems after converting to rain rate or ice mass. This can be 
accomplished by using Figs. 2 and 3 with a radar reflectivity - rain rate (Z-R) or other 
suitable relationship for a particular application. 

The reflectivity profiles verify that the precipitation index provides a measure of the 
precipitation and clouds in the vertical profile at nadir. When the API identifies neither 
rain nor cloud, there is usually no reflectivity detected by EDOP in Fig. 3a. Only in rare 
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cases does the reflectivity exceed 10 dBZ. Some non-precipitating clouds go undetected 
by API, but cloud identification is secondary to our goal of precipitation mapping. The 
cloud categories (1 and 2) also tend to have low (or sub-detectable) reflectivities (Fig. 3b). 
It is more common to find measurable reflectivity near the surface in category 2 (not 
shown), suggesting that shallow, very light rain is sometimes included. 

The shallow rain categories (3 to 5) do reliably indicate surface precipitation (Fig. 3c- 
d). These "shallow" rain profiles sometimes include an ice layer, but with low reflectivity 
values (i.e., small ice is unable to trigger the larger API categories). The higher 10.7 GHz 
emission thresholds in categories 4 and 5 result in larger low-level reflectivities from the 
liquid rain layer (Fig. 2). 

Categories 6 through 10, those having only the 85.5 GHz channel scattered, 
consistently include an ice layer detectable by EDOP (Fig. 3e-f). The reflectivities often 
decrease sharply from the liquid layer through the ice layer, indicating that any convection 
is weak. Radar bright bands are often present (except in category 10) and are indicative of 
stratiform rain. As intended, the liquid rain rates increase with increasing 10.7 GHz 
emission thresholds in categories 6 through 10 (Fig. 2). Reflectivities above the freezing 
level are similar for each of these categories, because they share the same 85.5 GHz 
scattering criteria. 

Categories 11 through 15 (i.e. those having the 37.1 GHz channel scattered) 
consistently have convective profiles with a deep layer of reflectivity greater than 20 dBZ 
(Fig. 2). Categories 16 to 18 (with the 19.35 GHz channel scattered) include the strongest 
convection, with 30+ dBZ radar echoes well above the freezing level (Fig. 2). In this last 
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set of categories (16-18), increasing 10.7 GHz TB no longer distinguish increasing rain 
rates. Instead, the lower 10.7 GHz thresholds for categories 16 and 17 sometimes 
distinguish greater scattering by large ice in these categories as compared to category 18. 
Categories 16 and 17 tend to have stronger convective profiles (greater reflectivity aloft) 
than category 18 (Fig. 2). This is an unintended result, but enables us to better resolve 
strong convection. The precipitation mass aloft in these categories attenuates the low 
level reflectivities. Even though attenuation has been accounted for following the alpha- 
adjustment technique (Iguchi and Meneghini 1994), comparison of liquid rain rates 
between categories 16 to 18 is not reliable. 

c. Merger of LIP information 

For many years, aircraft have routinely made measurements of electric fields 
associated with clouds, [e.g., Gunn and Parker. (1946); Blakeslee et al. (1989); Winn 
(1993)]. Some applications of aircraft electric field measurements have been made to the 
study of tropical cyclones [e.g., Orville et al. (1997), Simpson et al. (1998)]. For the ER- 
2, the vector electric field is derived using the outputs measured by the eight electric field 
mills installed on the aircraft The set of equations that relate these field mill outputs to 
the external electric field is represented as a matrix equation. Calibration of the field mill 
set on an aircraft involves the determination of the matrix coefficients using an iterative 
process (Mach and Koshak 2003). 

For this paper, calibrated electric field data are merged with the AMPR data by 
projecting the three-dimensional vector electric field onto the aircraft track as seen in 
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Figure 4. The three line plots in this figure represent the aircraft position projected onto 
the latitude/longitude, latitude/altitude, or longitude/altitude planes. The API data are 
mapped only on the latitude/longitude plane. The barbs on the aircraft tracks represent 
the two-dimensional projection of the vector electric field onto that plot at selected time 
intervals. The direction of the barb is the direction of the electric field (in that projection) 
while the length of the barb corresponds to the magnitude of the vector electric field (again 
in that projection). The resultant plot indicates the highest fields (where the barbs are the 
longest), as well as approximately where the sources of the fields are located (direction of 
the barbs). In general, for simple charge distributions, the electric field will point away 
(for positive charges) or towards (for negative charges) areas of charge. Note, while the 
low sample interval employed in this figure provides an excellent picture of the quasi- 
steady charge distributions in the clouds overflown, details of the transient field changes 
associated with lightning are not shown. 

4. Illustration - Hurricane Bonnie Case 

To demonstrate the relationships between the API, electric field, and radar 
reflectivity, a portion of the overflight of Hurricane Bonnie of 26 August 1998 is 
examined in detail. After becoming a hurricane on 22 August 1998 northeast of 
Hispaniola, Hurricane Bonnie made landfall near Wilmington, North Carolina on the 
afternoon of 26 August as a category 2 storm (Pasch et al. 2001). The NASA ER-2 
high altitude aircraft (carrying AMPR, EDOP, LIP, and other instruments) performed 
multiple overpasses of Hurricane Bonnie between 1 120 - 1720 UTC on 26 August. 
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Diverse precipitation structures were observed, allowing us to illustrate many API, 
electric field, and radar characteristics in a single example. As mentioned in the previous 
section, the API and electric field vectors are mapped in Fig. 4 between 1500-1600 UTC. 
This begins with a radial leg from the eye to the southwest quadrant, then a downwind leg 
from west to east, and an overflight of the eye with an exit to the northwest over land. 

The ER-2 flight track is overlaid on a 1501 UTC reflectivity image from the NOAA-42 P- 
3 lower fuselage radar in Fig. 5. The two figures show similar horizontal structure, 
considering that the precipitation field advects and evolves during the one-hour flight 
pattern. The greatest radar reflectivity values are seen in the western portion of the outer 
eyewall (not observed by this ER-2 flight segment but present in Fig. 5) and also in the 
southern and southeastern portions of die outer eyewall. The strongest electric fields and 
the largest API values (i.e., deepest and strongest convection) are observed while the ER- 
2 crosses the southeastern portion of the outer eyewall near 32.5N, 77.4 W (Fig. 4). In 
the southwestern portion of this flight segment, enhanced API and electric fields suggest 
strong convection to the left of the flight track. Several lightning flashes are detected, 
although they are not apparent with the data resolution plotted in Fig. 4. The suggested 
location of strong convection is consistent with the NOAA P-3 reflectivity patterns in 
the southern portion of Fig. 5, where enhanced eyewall reflectivity and some banding are 
seen just outside the eyewall. 

The southeast-to-northwest eyewall overpass in Fig. 4 is examined in further 
detail. Consider the vertical cross-section of reflectivity (Fig. 6), the vertical electric field 
and API (Fig. 7), nadir brightness temperatures (Fig. 8), and simulated radar reflectivity 
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(Fig. 9) derived from the median reflectivity profiles for each API value in Fig. 2. At far 
left (southeast) in Fig. 7, API = 1 suggests non-precipitating clouds. EDOP agrees, with 
only weak reflectivities (< 0 dBZ) around 6 and 12 km altitude. Following the flight 
track, next the API increases to values of 6 to 8, indicating rain with moderate ice. Indeed, 
the EDOP measures > 30 dBZ near the surface with > 10 dBZ reflectivity extending 
about 3 km above the bright band. The local maximum of API = 8 (heavy rain, moderate 
ice) coincides with the local reflectivity maximum (> 55 dBZ near the surface at x = 22 
km). The moderate ice categories are barely triggered, with the 85.5 GHz TB only a few 
degrees K less than the 37.1 GHz TB (Fig. 8). Because the ice scattering criteria are only 
minimally met, the reflectivity simulated by API in Fig. 9 overestimates the vertical 
extent of precipitation in this region. 

Continuing along the flight track, the API briefly decreases to 3 while reflectivity 
through the vertical column also decreases. The flight segment then encounters a thick 
anvil beginning around x=30 km, with echo tops reaching 16 km, API returning to values 
of 6 and 7, and the vertical electric field becoming slightly negative. This excursion of the 
electric field may be due to a weak positively charged layer near the top of the anvil, 
which is only ~4 km below the aircraft. 

API values between 6 and 9 vary with reflectivity in the rain layer. There is a 
close correspondence between API maxima (e.g., at x = 47 km; x = 75 km) and reflectivity 
maxima. The vertical electric field (E z ) becomes strongly positive and peaks at x = 90 km 
just before the eyewall reflectivity and scattering cores. This may be due to the sloping 
eyewall and the 20+ km altitude of the E z measurements. Peak E z at flight-level coincides 


70 


with peak reflectivity at -11-13 km altitude and the highest 10-20 dBZ echo tops). E* 
then decreases rapidly while the echo top heights also decrease. 

API increases to 14 and 15 (heavy rain with heavy ice) on the edge of the eyewall 
reflectivity core, and then increases to 16 through 18 (intense ice) over the core itself. 

This increase of API results from the lower frequency channels successively being 
scattered by larger graupel and/or hail. The peak 19.35 GHz scattering at x = 92 km is 
slightly offset from the peak 10.7 GHz emission at x = 95 km. This may be another 
result of the sloping eyewall, perhaps coupled with slight scattering in the 10.7 GHz 
channel at x = 92 km. The region with strongest electric field and most significant 
scattering also includes the strongest upper-level updrafts (Fig. 10). Future work will 
attempt to quantify this relationship. 

Inward (further right in the Figs. 6-10) from this outer eyewall is mostly shallow 
rain. Some of it is glaciated, but having low reflectivities above the bright band with 
minimal AMPR ice scattering signatures (values less than 6). AMPR does detect the 
inner eyewall on the southeast side at x = 150 km, with API = 7. Clouds are indicated 
inside the eye, with API values of 1 and 2 at x = 155-185 km. This is consistent with our 
visual observations from the NASA DC-8 aircraft, which was flying with the ER-2 during 
this mission. API suggests a broad region of shallow, light rain between x = 185-250 km 
and fails to detect the deep ice layer above the northwest inner eyewall at x = 200-220 
km. This feature is too weak for the API to handle properly; reflectivities are mostly 
below 20 dBZ both aloft and in the rain layer. There is a hint of 85.5 GHz scattering and 
some 37.1 GHz emission, but the algorithm requires more of either scattering or emission 
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in order to trigger an ice index (i.e., API 6 or greater, with TB85 < TB37). Consequently, 
the simulated reflectivity (Fig. 9) underestimates the vertical extent of this feature. 
Inclusion of a higher frequency channel (more sensitive to smaller ice) would likely help in 
situations such as this. 

Between x = 220-250 km, the API correctly identifies the shallow (~2 km) rain on 
the inner edge of the sloping outer-eyewall. The identification by API is qualitatively 
correct, but the simulated reflectivity overestimates the depth of die rain and 
underestimates the magnitude of the rain. EDOP vertical velocities show a shallow 
eye wall updraft (Fig. 10); the ice scattering design in the API is not particularly suited to 
detect this. On the far right (northwest side) of the cross-section, the rain increases in 
both depth and magnitude as API increases from 3 to 4 and then to 8. 

This cross-section was chosen because it demonstrates a wide variety of vertical 
structures, and includes all but the least common API values. Comparing API with 
reflectivity in this example, API does behave as qualitatively intended with very few 
exceptions. The simulated reflectivity shows structure similar to that measured by 
EDOP along the southeastern radial (the left half of the cross-section). In particular, the 
identification of localized rainfall maxima and the discrimination between different depths 
of precipitation is encouraging. The northwestern portion, particularly around x = 200- 
250 km, points out limitations of the API algorithm, as a deep layer of very light rain and 
ice is not distinguished from a shallow layer of heavier rain. 


7.7 


5. Summary 


An oceanic precipitation screening technique that combines the information 
content of the four AMPR frequencies at a given data pixel into one precipitation index is 
presented. The technique, which has been verified with EDOP and NOAA P-3 data, 
shows promise as a computationally easy rainfall mapping tool suitable for application to 
high spatial and temporal resolution airborne data. Merger of die precipitation index with 
three-dimensional electric field data readily identifies the convective strength of embedded 
cells within precipitation systems. Further study using this type of analysis will examine 
the other CAMEX precipitation cases to quantify the relationship between lightning and 
microwave information as a surrogate indicator of convective strength. A more extensive 
examination of the NOAA P-3 radar information and the CAMEX microphysical data 
will be conducted to explore the feasibility of adding a rain rate conversion algorithm to 
the API screening technique for use as quantitative precipitation estimation tool. 

The synergy of the AMPR, EDOP, and LIP data sets has been presented here not 
only as research tool for those interested in hurricane studies or as a validation tool for 
those developing satellite rainfall algorithms, but also as an example of how airborne 
information may be merged into real-time observational products. Future concepts for 
Earth observation include adding airborne platforms such as uninhabited aerial vehicles or 
ultra-long duration balloons into a mixture of spacebome and surface-based assets 
comprising a flexible, adaptive global observation network. Within these types of 
frameworks, an airborne vehicle could be positioned to provide high spatial and temporal 
coverage of a critical weather event in concert with spacebome and surface 
instrumentation so that the best combination of information is used for observation and 
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prediction of the event outcome. As technical development of airborne platforms for this 
type of use progresses, appropriate airborne instrumentation and data algorithms should 
be identified that provide the maximum amount of information using the most feasible 
airborne payload for a given application. This study presents instrument candidates that 
could be used for high altitude monitoring of precipitation type and convective strength 
for tropical cyclone and other precipitation systems. The research and operational 
communities should also examine many other types of instruments and flight altitudes in 
order to choose the optimal mixture of observations. 
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Table 1 The AMPR precipitation index descriptions, criteria, and representative rain 
rates. Rain rates were estimated by applying the Jorgensen and Willis (1982) 
Z=300R’ 35 relationship to the median 1 km altitude reflectivity value for each 
FDOP/AMPR vertical nrnfile 


Index 

Color 

Description 

Criteria 

Dn* 

Rain 

rate 

0 

a 

Clear 

TB10 < 160 and TB37 < 215 

n/a 

0 

1 

n 

Light-to-Moderate Cloud 

TB19 > 190 or TB85 > 260 

n/a 

0 

2 

63 

Moderate-to-Heavy Cloud 

TB85 > 270 

n/a 

0 

& 

II 

f 

Ice Level 0, indices 0-5 

TB85 > TB37 or TB85 > 275 and 


3 


Level 1 Rain 

TB10 > 160 or TB37 > 215 

.40 

0.4 

4 


Level 2 Rain 

TB10 > 175 

.005 

5 

5 


Level 3 Rain 

TB10 > 200 

.001 

7 


Ice=lE3 

Ice Level 1, indices 6-10 

TB85 < TB37 and TB85 < 275 and 


6 


Level 1 Rain 

TB10 > 160 or TB37 >215 

.27 

2 

7 


Level 2 Rain 

TB10 > 175 

.14 

5 

8 


Level 3 Rain 

TB10 > 200 

.05 

8 

9 


Level 4 Rain 

...... . 

TB10 > 225 

.03 

13 

10 


Level 5 Rain 

TB10 > 250 

.02 

19 

□ 


Ice Level 2, indices 11 - 15 

Ice 1 and TB37 < TB19 and TB37 < 260 and 


11 

■ 

Level 1 Rain 

TB10 > 160 or TB37 > 215 

.002 

3 

12 

■ 

Level 2 Rain 

TB10 > 175 

.02 

4 

13 

m 

Level 3 Rain 

TB10 > 200 

.03 

7 

14 


Level 4 Rain 

TB10 > 225 

.02 

13** 

15 

n 

Level 5 Rain 

TB10 > 250 

.01 

19** 

■ 

Ice=m 

Ice Level 3, indices 16 - 18 

Ice 2 and TB19 < TB10 and 


16 


Level 4 Rain 

TB10 > 225 

<.001 

NA 

17 


Level 5 Rain 

TB10 > 250 

.005 

27** 

18 


Level 6 Rain 

TB10 > 275 

.003 

37** 


is the fraction of the index to the total of all rain-only samples (t| =231042). 

** estimated at 4 km altitude due to attenuation uncertainty 

Notes: 

Ice Level 0 (No Ice) and Levels 4-6 Rain do not significantly occur and are included in Index 5. 

Ice Level 1 (Moderate Ice) and Level 6 Rain does not significantly occur and is included in Index 10. 
Ice Level 2 (Heavy Ice) and Level 6 Rain does not significantly occur and is included in Index 15. 

Ice Level 3 (Intense Ice) and Levels 1-3 Rain do not significantly occur and are not included. 




Brightness Temperature in Kelvin 


Fig. 1. AMPR brightness temperatures at 10.7, 19.35, 37.1, and 85.5 GHz for 
Hurricane Bonnie at 1500-1600 UTC on 26 August 1998. Low brightness 
temperatures are magenta and blue, high brightness temperatures are yellow and red. 
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Fig. 4. Horizontal mapping of API and projections of 3-D electric field and aircraft 
location onto (A) longitude-altitude plane, (B) latitude-altitude plane, and (C) 
longitude-latitude plane. The "s’ denotes start of aircraft track (1500 UTC) and the 'e' 
denotes end of aircraft track (1600 UTC). API color scale as in Table 1. Electric field 
and aircraft location are plotted as if projected onto three sides of a box. The 
projections of the electric field onto the (A) longitude-altitude, (B) latitude-altitude, 
and (C) latitude-longitude planes are plotted as barbs originating at the aircraft 
location. A barb extending 1 km above the aircraft track denotes a +1 kV m' 1 (positive 
charge below the aircraft) electric field. Note that the vertical component dominates 
the longitudinal and latitudinal components of electric field in (A) and (B). Barb 
lengths are scaled by a factor of 10 in (C) because the horizontal components of 
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Fig. 5. Reflectivity scan from NOAA-42 lower fuselage radar in eye of Hurricane 
Bonnie at 1500 UTC 26 August 1998. ER-2 flight track between 1500-1600 UTC is 
overlaid. 
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Fig. 6. Vertical cross-section of nadir EDOP reflectivity, 1535 - 1557 UTC 26 August 
1998. Cross section extends from 170 km southeast of the center of Hurricane Bonnie 
(left) to 1 10 km northwest of the center (right). Contours every 5 dBZ; colors as in Fig. 
9. 


36 



Hor izo nta l dfatance (Km) 


Fig. 7. Nadir API (blue) and vertical component of electric field (red) coincident with 
Fig. 6. 
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Fig. 8. Nadir AMPR brightness temperatures at 10.7, 19.35, 37.1, 85.5 GHz 
coincident with Fig. 6. 
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Fig. 9. Radar reflectivity simulated by API, coincident with Fig. 6. Simulated 
reflectivity is the convolution of observed API in Fig. 7 and the median reflectivity 
profiles in Fig. 2. Black contours are measured EDOP reflectivity from Fig. 6. 




Fig. 10. Vertical cross section of nadir EDOP vertical velocity 
coincident with Fig. 6. Flydrometeor fallspeeds removed following 
Heymsfield et al. (1999). 
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LIST OF CAPTIONS 


Table 1 The AMPR precipitation index descriptions, criteria, and representative rain 
rates. Rain rates were estimated by applying the Jorgensen and Willis (1982) Z=300R 1 35 
relationship to the median 1 km altitude reflectivity value for each EDOP/AMPR vertical 
profile. 

Fig. 1. AMPR brightness temperatures at 10.7, 19.35, 37.1, and 85.5 GHz for Hurricane 
Bonnie at 1500-1600 UTC on 26 August 1998. Low brightness temperatures are magenta 
and blue, high brightness temperatures are yellow and red. 

Fig. 2. Vertical profiles of median EDOP reflectivity for all nadir pixels, sorted by API. 
API=1 is omitted because it is subzero at all heights. API=16 is omitted due to 
insufficient sample size. 

Fig. 3. Cumulative Density Function of EDOP reflectivity as a function of height for 
selected API values. 

Fig. 4. Horizontal mapping of API and projections of 3-D electric field and aircraft 
location onto (A) longitude-altitude plane, (B) latitude-altitude plane, and (C) longitude- 
latitude plane. The ‘s’ denotes start of aircraft track (1500 UTC) and the ‘e’ denotes end 
of aircraft track (1600 UTC). API color scale as in Table 1. Electric field and aircraft 
location are plotted as if projected onto three sides of a box. The projections of the 
electric field onto the (A) longitude-altitude, (B) latitude-altitude, and (C) latitude- 
longitude planes are plotted as barbs originating at the aircraft location. A barb extending 1 
km above the aircraft track denotes a +1 kV m" 1 (positive charge below the aircraft) 
electric field. Note that the vertical component dominates the longitudinal and latitudinal 
components of electric field in (A) and (B). Barb lengths are scaled by a factor of 10 in 
(C) because the horizontal components of electric field are so small. 

Fig. 5. Reflectivity scan from NOAA-42 lower fuselage radar in eye of Hurricane Bonnie 
at 1500 UTC 26 August 1998. ER-2 flight track between 1500-1600 UTC is overlaid. 

Fig. 6. Vertical cross-section of nadir EDOP reflectivity, 1535 - 1557 UTC 26 August 
1998. Cross section extends from 170 km southeast of the center of Hurricane Bonnie 
(left) to 1 10 km northwest of the center (right). Contours every 5 dBZ; colors as in Fig. 
9. 

Fig. 7. Nadir API (blue) and vertical component of electric field (red) coincident with Fig. 

6 . 
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Fig. 8. Nadir AMPR brightness temperatures at 10.7, 19.35, 37.1, 85.5 GHz coincident 
with Fig. 6. 

Fig. 9. Radar reflectivity simulated by API, coincident with Fig. 6. Simulated reflectivity 
is the convolution of observed API in Fig. 7 and the median reflectivity profiles in Fig. 2. 
Black contours are measured EDOP reflectivity from Fig. 6. 

Fig. 10. Vertical cross section of nadir EDOP vertical velocity 
coincident with Fig. 6. Hydrometeor fall speeds removed following 
Heymsfield et al. (1999). 
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